Ballistic spin-polarized transport and Rashba spin precession in semiconductor 

nanowires 
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We present numerical calculations of the ballistic spin-transport properties of quasi-one- 
dimensional wires in the presence of the spin-orbit (Rashba) interaction. A tight-binding analog 
of the Rashba Hamiltonian which models the Rashba effect is used. By varying the robustness of 
the Rashba coupling and the width of the wire, weak and strong coupling regimes are identified. Per- 
fect electron spin-modulation is found for the former regime, regardless of the incident Fermi energy 
and mode number. In the latter however, the spin-conductance has a strong energy dependence due 
to a nontrivial subband intermixing induced by the strong Rashba coupling. This would imply a 
strong suppression of the spin-modulation at higher temperatures and source-drain voltages. The 
results may be of relevance for the implementation of quasi-one-dimensional spin transistor devices. 



PACS numbers: 72.25.Dc, 73.23.Ad, 73.63.Nm 
I. INTRODUCTION 

The electron spin precession phenomena at zero mag- 
netic field induced by a variable spin-orbit interaction 
in quasi-two-dimensional electron gas (2DEG) systems, 
was first elucidated by Datta and DaaH as the basic 
principle for the realization of a novel electronic de- 
vice, the spin-transistor. The underlying idea is to 
drive a modulated spin-polarized current (spin-inject and 
detect) entirely electrically, combining just ferromag- 
netic metals and semiconductor materials. For this, the 
spin-precession ia-jcpntrolled via the spin-orbit coupling 
(Rashba-couplingtffl) associated with the interfacial elec- 
tric fields present in the quantum well that contains the 
2DEG channel. The tuning of the Rashba coupling by 
an external gate voltage, was recently demonstrated in 
different semiconductors by Nitta et ala and othersoQ, 
and more recently by Grundlem applying a back gate 
voltage while the carrier density was kept constant. It 
has been also achieved in a p-type InAs semiconductor 
by Matsuyama et alB 

Although spin injection has already been reported fram. 
ferromagnetic metals to InAs based semiconductors J13'EJ 
the spin polarization signatures reported are about 1% or 
less, making the results very controversial. Such low effi- 
ciency can be also attributed to extraneous effects, such 
as the local Hall effect ,E30 and to resistance mismatch 
at the interfaces between the ferromagnetic metal and 
the semiconductor .113 Very recently, one of us has pro- 
posed that the latter problem can be circumvented by 
growing atomically ordered and appropriately oriented 
interfaces of ferromagnetic metals and suitable semicon- 
ductors, that act as perfect spin-filters and make injection 
of up to 100% spin-polarized electrons into semiconduc- 
tors possible in principle. Ej If this prediction is confirmed 
experimentally, a major obstacle to the spin injection ex- 
periments will be overcome. 



On the other hand, another crucial prerequisite to 
having an overall strong spin current modulation, is 
to restrict the angular distribution of electrons in the 
2DEG by imposing a strong enough, transverse confin- 
ing potentials This was the original proposal of Datta 
and Das. It was argued that if W is the width of 
the transverse confining potential well, the condition 
W <C h 2 /am*, should be satisfied for the intersubband 
mixing to be negligible, which would ensure a perfect 
spin-current modulation. If we choose a typical channel 
in an Ino.53Gao.47As semiconductor, with m* = 0.042 m— 
and spin-orbit coupling constant a = 1.5 x 10 _11 eVm,tl 
this implies W <C 0.38/Ltm. This width is of the order 
or much smaller than the characteristics lateral widths 
of the 2DEG channels patterned in the current exper- 
imental studies of spin-injection.EJiij This would sug- 
gest that in order to achieve better results in the spin- 
injection modulation, the best choice would be ballistic 
quasi-one dimensional systems (by introducing split gates 
for instance), where just a few populated subbands are 
allowed, rather than the quasi-two-dimensional situation, 
in which many propagation channels exist (Hall-bar ex- 
periments), with a concomitant small subband spacing 
comparable to the zero field spin-splitting energy of the 
2DEG. 

Recently, Moroz and BarnesQ in a theoretical study 
of the effect of the spin-orbit interaction on the bal- 
listic conductance and the subband structure of quasi- 
one-dimensional electron systems, showed that a dras- 
tic change in the fc-dependence of the subband spectrum 
occurs with respect to the purely 2DEG system when 
relatively strong spin-orbit coupling is considered. This 
yields additional subband extrema and subband anti- 
crossings, as well as anomalous peaks in the conductance 
of the Q1DEG. To what extent these effects can influence 
the behavior of a quasi-one-dimensional spin-modulator 
device has not been investigated, and this is the aim of 
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the present work. 

In this paper we investigate the effect of the strength 
of the Rashba spin-orbit coupling on the spin transport 
properties of narrow quantum wires. We find it conve- 
nient for this purpose to work in a simple tight-binding 
approach in which an homologous version of the Rashba 
spin-orbit coupling is employed. In particular, we find 
that the spin-orbit interaction induces dramatic quali- 
tative changes in the spin-polarized current transmitted 
through Q1DEG systems, provided that a strong spin- 
orbit coupling is present. A strong dependence of the 
spin-conductance on the incident Fermi energy is found 
to be correlated with subband mixing induced by a robust 
spin-orbit coupling. This dependence can significantly 
suppress the spin-modulation at finite temperatures and 
or bias voltages. These results may have important impli- 
cations for prospective quasi-one-dimensional spin injec- 
tor devices at room temperature or under large applied 
voltages. 

The remainder of the article is organized as follows: 
In Sec. II the theoretical approach is developed, start- 
ing with a brief summary of the relevant features of the 
Rashba Hamiltonian and its induced spin precession ef- 
fect. A tight-binding model for the Rashba Hamiltonian 
is also presented here. A brief description of the con- 
ductance calculation is given at the end of the section. 
The numerical results and main conclusions are given in 
Sec. Ill, and finally the criteria that distinguish between 
the weak and strong spin-orbit coupling regimes and the 
method used to obtain the subband spectrum are out- 
lined in the Appendices. 

II. THEORETICAL MODEL 

A. Hamiltonian for the Rashba effect 

In the absence of a magnetic field, the spin degeneracy 
of the 2DEG energy bands at k ^ is lifted by the cou- 
pling of the electron spin with its orbital motion. This 
coupling arises because of the inversion asymmetry of the 
potential that confines the 2DEG system. The spin-split 
dispersion involves a linear term io-ik, as was first in- 
troduced by Bychkov and Rashba.Em The mechanism is 
popularly referred to as the Rashba-effect. The spin-orbit 
(Rashba) model is described by the Hamiltonian 

H so = ^(<? x p) z = ia(a y -^- - a x^) (1) 

Here the z axis is chosen perpendicular to the 2DEG 
system (lying in the x — y plane), a is the spin-orbit 
coupling constant (Rashba parameter) which is sample 
dependent, and is proportional to the interface electric 
field, a — (<j x , <7 y , tj z ) denotes the spin Pauli matri- 
ces, and p is the momentum operator. The experimen- 
tal values of a for different materials range from about 
6xl0~ 12 eVmat electron densities of n — 0.7xl0 12 cm -2 



to 3.0 x kQj-JieVm at electron densities of n = 2 x 
10 12 cm" 2 ™ 

The Rashba Hamiltonian CL), which is derivable from 
group theoretical argumentst3, is invariant under time 
reversal, that is [T, H so ] = 0. The time reversal operator 
is represented here by T — ia y C, with C the complex 
conjugation operator. Since the degeneracy of the elec- 
tronic states at k = can be only lifted if the time re- 
versal symmetry of the system is broken, it follows that 
the Rashba Hamiltonian, (due to its time reversal invari- 
ance), can not produce a spontaneous spin polarization 
of the electron states. Nevertheless, as mentioned earlier, 
it is capable of removing the spin degeneracy for k / 0. 
This is made clear by noticing that the total effective 
mass Hamiltonian for a 2DEG system, as a result of the 
Rashba effect [Eq. (1)], has the form 

H= ( *&( k Z + k D iak x + ak y \ 
\ ~iak x + ak y ^(fc 2 + fc 2 ) J ' 

where H = H +H so , with H a being the electronic kinetic 
energy part in the absence of the Rashba effect. Clearly, 
the Hamiltonian H produces two separate branches for 
the electron states, 

h 2 

E(k) = k 2 ±ak, (3) 

here k — |k| is the magnitude of the 2D wavevector in 
the 2DEG plane. Since the spin-orbit coupling a de- 
pends on the interface electric field, it is possible to tune 
the strength of the splitting between the two branches by 
applying an external gate voltage, which will alters the 
net effective electric field at thfeijiterface, as has been 
verified in several experimentsJStm 

B. The Rashba spin-precession effect 

The Rashba effect is the basis for the proposed, and 
yet to hue implementated, Datta and Das spin-modulator 
device.clln this device, a spin-polarized current is injected 
from a ferromagnetic material into a 2DEG at a inver- 
sion layer (formed at a semiconductor heterojunction) 
and then collected by a second ferromagnetic material 
(Fig. 1(a)). In basic terms, the idea is that the Rashba 
effect will induce a spin-precession of the electrons mov- 
ing parallel to the interface, rotating them with respect 
to the magnetization direction of the second ferromag- 
net (collector). Then by adding a gate voltage the net 
effective electric field (and hence, the spin-orbit interac- 
tion) at the interface can be modified, tuning the spin 
precession, and therefore, the transmitted spin-polarized 
current is modulated accordingly.EI 

We mentioned in the introduction that an important 
prerequisite to having an overall strong spin current mod- 
ulation that the angular distribution of the 2DEG be 
restricted by imposing a transverse confining potential. 
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Bearing this in mind, we will now summarize the essential 
physics of the spin-precession effect in a Q1DEG system 
for the case of weak spin-orbit coupling. The summary 
will define the basic conceptual framework that will be 
needed to understand our results in the general multi- 
channel case with arbitrary spin-orbit coupling strength 
and will also establish the notation to be used in the 
remainder of the paper. 

Consider a Q1DEG system which is defined by ap- 
plying split gates to a 2DEG in a semiconductor het- 
erostructure (In a; Gai_ 3 ;As/InyAli_yAs for instance, see 
Fig. 1(b)). Due to the confining potential V(x), the 
electron motion will be quantized in the x direction, 
Fig. 1(c). Following Datta and DasEI, let's assume 
that the Rashba spin-orbit interaction H so is sufficiently 
weak that its effect can be incorporated pertubatively. 
For such a case, the unperturbed (a = 0) Hamiltonian 
will satisfy H \n, a) — E°\n, a), where the eigenvalues 
are given by E° = E n + h 2 k 2 /2m*, with n denoting 
the subband index. The unperturbed spin degenerate 
eigenstates have the form, \n,a) — ► e lhyV (f) n (x)\a) , where 
a =1,1 with the definitions of the spinors | f) = ( □ ), 
and | 4) = ( i )• Note that the <fi n (x) are the solutions of 



2m* dx 2 



+ V(x) (j) n {x) = E n (j) n (x). 



(4) 



We seek the eigenvalues for the system in which we 
have a non-zero and weak spin-orbit interaction, a ^ 
(i.e. H so ^ 0). That is for H = H a + H so . From degen- 
erate perturbation theory, we obtain (to zeroth oder) the 
following system of equations for each subband n, 



( E n + {HsoYnn ~ E )a° na + ) (H so 



aa a° , 

nn ^na' 



(5) 



with a° na the zeroth order expansion for the coefficients 
a nrT used to expand the perturbed states in terms of the 
known unperturbed states |n, a). The result in Eq. (5) 
is valid as long the condition 



(H ao ) 



1 1 in 



(E° m - E° n ) 



< 1, 



(6) 



\crcr 
I nm 



for m ^ ii has been fulfilled, where {H so ) c r 
(n, <j\H so \m, a 1 ) are the matrix elements that intermix 
the different subbands and spin-states in the perturbed 
system. Explicitly they are given by 



{H so ) 



Tl 

nm 



akydnm + a(n\ — \m) 

{H so )nm = akySnrn ~ J-H m } 



{Hso) n ln 



(H 



a 

so Jnrn 



o. 



(7) 

(8) 
(9) 



Clearly the perturbation is non-diagonal in the spin 
index and it is linear in k y . If we suppose that the trans- 
verse confining potential V(x) has reflection symmetry 
in x, eqs. (7) and (8) will reduce to 



{H eo )H = (H so [ 



IT 

nn 



ak 1t 



(10) 



Inserting the results from (9) and (10) into equation (5) 
we get for each channel n, 



E° n -E ak y \ < T) 
ak y E°-Ej{a° nl 1 



which yields the eigenvalues 



E ± (k y ) = E n + —k 2 ± aky. 

2m* • 



(11) 



(12) 



Eq. (12) shows that for Q1DEG systems, the Rashba 
spin-orbit interaction introduces (to zeroth order) , a lift- 
ing of the spin degeneracy for each subband state n. The 
nature of the splitting is such that it allows electrons 
with the same energy to have different wave vectors, (k y i 
and k V 2), that is E + (k y i) = E~(k y 2), where k y \ is the 
wave vector associated with the subband E + with eigen- 
vector ( J ) , whereas k y i represents the wave vector as- 
sociated with the subband E~ with eigenvector f \ ) . 
Therefore, if we were to drive spin-up polarized electrons 
into the Q1DEG using a ferromagnetic material as in 
the Datta and Das device (Fig. 1(a)), the wave emerg- 
ing from the semiconductor wire would be represented 
as ijj =( \ \ e lkylL + ( _\ ~)e lky2L , where L is the length 
of the semiconductor wire. Therefore the probability of 
detecting a spin up ( J ) electron at the collector contact 
would be proportional to 



|((10)|^)| 2 =4cos i 



(ky 2 - kyl)L 



(13) 



whereas if the collector contact is magnetized such that 
it detects only spin-down polarized electrons, the proba- 
bility of detecting a spin-down ( ° ) electron will be pro- 
portional to 



|<(0lM| 2 = 4sim 



(ky2 ~ k y l)L 



(14) 



The results (13) and (14) are very important; they im- 
ply that if we inject and collect spin-polarized (up or 
down) electrons into (and from) a Q1DEG system, the 
Rashba-effect will produce a modulation of the trans- 
mitted current at drain contact with a differential phase 
shift given by A6 — Ak y L, where Ak y — k y 2 — k y i . 
In other words, the Rashba effect would induce a spin 
precession of the transmitted electrons with a phase shift 
A9 with respect to those injected at the ferromagnetic 
emitter. From Eq. (12) is straightforward to determine 
Aky. Since iS+^i) = E~(k y2 ) we have that 

h 2 

E+(k yl ) - E (kyi) = ^—Mi - K 2 ) + a(kyi + k y2 ) 
= (15) 
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which leads to k y 2 — k y i = 2m*a/K 2 . Therefore the dif- 
ferential phase shift A9 can be written aso, 



» /. 2m* 
A9 = -^-aL, 



(16) 



which is proportional to strength of the Rashba param- 
eter a and to the separation L between the magnetic 
contacts or the length of the Q1DEG system. Then by 
applying a back gate voltage bias, the Rashba parameter 
can be varied in principle, and hence, the degree of the 
electron spin precession would be tuned correspondingly. 

Eq. (16) suggests that spin-current modulation can 
be attained in semiconductor nanowires, regardless of 
the mode number and large applied bias, as indeed 
Datta and Das concluded.EJ Notice as well that the phase 
shift is energy independent. We should recall here that 
these results were obtained with the premise that the 
Rashba spin-orbit interaction was weak enough, such 
that the confinement energy was much larger than the 
spin- splitting energy induced by the Rashba effect, and 
therefore, the intersubband mixing was neglected. In 
other words, whenever the condition (6) is satisfied. For 
the hard wall confining potential considered here this 
means 



(H so ) 

III 



am*W 



< 1 



(17) 



which would imply W <C h 2 /am* and gives us a rough 
upper limit of the width of the confinement for the inter- 
subband mixing to be neglected. Typical values for m* 
and a yield widths of W <C 0.38 /im. We will arrive at 
the similar conclusions using parabolic potential for the 
transverse confinement V(x). 

The criterion (17) for the validity of the pertubative 
theory discussed above is quite restrictive since it requires 
both weak spin-orbit coupling and a narrow wire. Thus in 
general the intersubband mixing that has been neglected 
in the above discussion can be important, and the sim- 
ple energy-independent expression (16) for the differen- 
tial phase shift may then no longer apply. This results in 
important qualitative changes in the electron spin preces- 
sion and in the behavior of the spin-modulation devices 
as will we discuss below. 

In the following section we present a tight-binding ana- 
log of the Rashba Hamiltonian that can be used to study 
spin transport in cases where intersubband coupling is 
important and that pertubative theory fails, and also in 
more general geometries which do not lead themselves 
readily to analytic solution. We then proceed to solve the 
multi-channel scattering and spin-dependent transport 
problem exactly (numerically) within the tight-binding 
formalism in a variety of situations where intersubband 
mixing is significant and needs to be considered. 



C. Rashba Hamiltonian: A tight binding model 

In this section, the spin-orbit interaction given by the 
Rashba Hamiltonian is reformulated within the tight- 
binding approach in a lattice model. We consider a quasi- 
one dimensional wire, which is assumed to be infinitely 
long in the propagation direction. The wire is represented 
by a two dimensional grid with lattice constant a. We 
choose the coordinate system such that the x axis, with 
N x lattice sites, is in the transverse direction, while the y 
axis, with N y lattice sites is in the longitudinal direction, 
Fig. 2(a). 

We assume only nearest neighbor spin-dependent in- 
teractions for the Rashba perturbation. For the purpose 
of the present model, the localized site orbitals will be 
assumed to have the symmetry of s-states. Then the 
tight-binding analog of eq. (1) takes the form 

Hlo = - t so ^2 X/ (pl+l,m,(r'^' <T v)<r,<r'^'t,m,a 
a,a f l,m 

~ Cl^yiia^Ct^) + H.c. (18) 

with an isotropic nearest-neighbor transfer integral t so 
which measures the strength of the Rashba spin-orbit 
interaction and will be shown below to have the value 
t S o = a/2a, C\ represents the electron creation oper- 
ator at site (l,m) with spin state a, (a =T, !)• Here a x 
and <j y are the Pauli spin matrices. Notice that t so > 0, 
and in principle it can be tuned by an external electric 
field. 

This Hamiltonian is formally -similar to that studied 
previously by Ando and TamuraEEl in the context of con- 
ductance fluctuations and localization in quantum wires 
with boundary-roughness scattering. We emphasize that 
in the present model no roughness scattering is present, 
and that here the physical origin of spin-orbit scattering 
is the asymmetry of the electric field in the quantum well 
that contains the 2DEG. 

We divide the wire into three main regions. In two of 
these (I and III in Fig. 2(a)) which are near the ferromag- 
netic source and drain, the spin-orbit hopping parameter 
t so is set to zero. We assume that the semiconductor 
interface at which the Rashba effect occurs does not ex- 
tend into these two regions and that t so = there for 
that reason. In the middle region (II) the spin-orbit cou- 
pling is finite, (t so ^ 0) at the semiconductor interface. 
In the actual calculation this region (II) is further di- 
vided into three, to introduce the spin-orbit interaction 
adiabatically. Therefore, the full Hamiltonian reads 



tb 



H = H + V + H 
where the spin diagonal parts of H are given by 



(19) 
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and 



£,m,cr 



(20) 



(21) 



with e the on-site energy, t is the hopping energy (t = 
~h 2 /2m*a 2 ), and V/ )7 „ is an additional confining poten- 
tial. The full Hamiltonian H will have eigenvectors given 
by 



(22) 



where we have defined the spinors, 

= xl, m ( J ) . and = 4 m ( J ) (23) 

with xf m = (Zmal^), xfm being the electronic wave 
function at site (l,m) and in spin state a =\, J.. In (22), 
\l m a) = Cj ma \0), with |0) denoting the vacuum, and 
we assume (V ml a'\l m a) = 5i,i>& m ,m'5a,o' 



We establish the correspondence between Eqs. (18) 
and (1) and determine the value of t so by noticing that 
for a 2DEG system with plane wave solutions of the 
form e i( k * al + k v am ) ; the Hamiltonian H yields the two- 
dimensional tight-binding eigenvalues, 



E(k\\) = E (k\\) ± 2t so ^ sin 2 (k x a) + sin 2 (k y a) (24) 

with ^(kn) — e + 2t[cos(k x a) + cos(k y a)] the tight- 
binding conventional subband dispersion for a 2DEG, 
with k|| = (k x ,k y ,0). For small k x a and k y a (setting 
the on-site energy e equal to —it), Eq. (24) reduces to 



E 



2m* 



{k 2 x + k 2 y )± 2at so Jk[+ 



k y 



(25) 



Note that Eq. (25) is just the continuum subband dis- 
persion (3), and will correspond to the expected Rashbi 
subband linear splitting with the definition a = 2at so H. 



D. Spin transport calculation 

For the study of the spin-dependent transport in a 
Q1DEG system, the physical model we bear in mind is 
as follows: Consider a Q1DEG spin-modulator device as 
shown in Fig. 2(b). The device has two independent 
ferromagnetic source contacts with magnetizations ori- 
ented such that one of them can emit only spin-up polar- 
ized electrons (contact SF^)), whereas the second source 
contact SF^-*, can emit only spin-down polarized elec- 
trons. Likewise, two independent ferromagnetic drain 



contacts DF(T) and DF^ are attached at the opposite 
end of the device which are able to detect just spin-up 
and spin-down polarized electrons, respectively. A per- 
fect ohmic contact between the ferromagnetic materials 
and the semiconductor is assumed. A back gate under- 
neath the device and directly below the Q1DEG channel 
will control the spin precession (through the Rashba ef- 
fect) of the injected electrons. We will suppose that the 
device is set up such that spin polarized electrons are 
launched either from the SF^) or SF^) contacts, while 
both spin orientations (f and |) are being absorbed at 
the drain contacts, ensuring in this way a spin-resolved 
measurement. In addition it is supposed that the elec- 
tron spin reorientation due to defect scattering can be 
neglected, i.e., the spin relaxation time is much longer 
than the electron dwell time in the device. 

We now discuss the approach used in this work 
to study spin-transport in the Q1DEG system de- 
scribed above, specifically, for the calculation of the 
spin-conductance. The spin-transport problem was 
solved numerically through the use of the spin-dependent 
Lippman-Schwinger equation, 



|*) = \$)+G (E)U\V) 



(26) 



where |<&) is the unperturbed wave function, i.e. an 
eigenstate for H , while G (E) = (E + ie — i? ) _1 is 
the Green's function for the system in the absence of any 
kind of scattering. Here U = V + H so , represents the 
scattering part, with V the scattering due to the con- 
finement potential, and H so the spin-dependent part due 
to the Rashba coupling in the semiconductor wire. The 
unperturbed functions, |$) and G (E) are, known ana- 
lytically, see for instance Nonoyama et alr3 Notice that 
the unperturbed Green's function is also diagonal in spin 
index, that is 



(l q m q cr q \G \lma) = (l q m q \G \lm)8 atTq , 



(27) 



therefore, for a wave function at any site (l q ,m q ) and 
state cr q , the Lippman-Schwinger equation for this sys- 
tem can be rewritten as follows 



= $ £m,+ E [{l q m q \G \l m)V lm ^l m S ataq 

- t S0 (l q m q \G \l + 1 m)(ia y ) atrTq ^l m 
+ t S0 (l q m q \G \l m + l){ia x ) 

- t S0 (l q m q \G \l m)(ia y )* a ^>f +1 



(28) 



+ t S0 (l q m q \G \l m)(ia x )*^f m+1 



By solving the coupled linear equations resulting from 
the equation above, it is possible to determine the 
perturbed wave function associated with the complete 
Hamiltonian H at any site of interest. 

Once the wave functions for each state a are known 
at the ferromagnetic contacts (regions I and III in Fig. 
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2(a)) the spin-dependent .conductance is obtained within 
the Landauer frameworkEj, which is appropriate if the 
electron-electron interactions are unimportant £3 In our 
calculations, only spin-up polarized electrons are injected 
from the emitter (region I in Fig. 2(a)) into the spin- 
orbit interaction region (where the spin precession of the 
incident electron is induced). At the collector region, the 
transmitted current is calculated separately for each spin, 
modeling the pair of drain contacts DF^' and DF^ in 
Fig. 2(b). Therefore, in general we will have two contri- 
butions for the net transmitted current, one coming from 
spin-up electrons that arrive at the collector (region III) , 
and the other one from the collected spin-down electrons. 
We therefore define the conductances G^ and G^ associ- 
ated with the currents flowing between SF^) and DF^- 1 , 
and DFw, respectively, by 

V 

and 

V 

It is important to note that the spin-down partial con- 
ductance at the collector ferromagnet G^ arises due to the 
induced spin precession of the incident spin-up polarized 
electrons and since no spin-down polarized electrons are 
being injected. In Eq. (29) is the partial transmission 
probability (summed over the incident channels) that an 
incident electron (injected from SF^- 1 ) with spin a =j, 
is transmitted to the right ferromagnet DF^) ; v denotes 
the outgoing transmitted mode, while is the partial 
transmission probability that an incident electron with 
the same spin a =f, is transmitted as a spin-down elec- 
tron, and measured at DF^ . The partial transmission 
probability is given explicitly by 

< = £(§)i^ M 

where v% and v° are the outgoing and incident electron 
velocities at the Fermi energy with spin a and modes v 
and /i of the drain and source contacts, respectively. In 
(31) tj is the partial transmission amplitude that an in- 
cident electron with spin a —\ and mode fj, is transmitted 
in the v mode with the same spin state a =f, and t+„ is 
the partial transmission amplitude that an incident elec- 
tron with spin a =\ and mode \i is transmitted in the v 
mode with the opposite spin state, i.e. a =J, . 

III. NUMERICAL RESULTS AND THEIR 
IMPLICATIONS 

In order to study the dependence of the spin- 
conductance on the strength of the spin-orbit interac- 
tion, we will consider two cases, namely, weak and strong 



coupling. The criteria that distinguish these two cases 
are obtained in the appendix A; the physical considera- 
tions are essentially as follows: Since in the multichan- 
nel scattering process the eigenstates of the full Hamilto- 
nian are (in general) linear combinations of the different 
spin-subbands (due to the Rashba-term) , therefore, in 
a pertubative sense, the contribution of the mixing of 
the spin-subbands should be negligible as long the sub- 
band spacing (separation in energy) AEw = E m — E n , 
is much greater than the subband intermixing energy, 
Eq. (6). However, if the confinement energy and/or the 
spin-orbit coupling are of the same order as the energy 
shift introduced by the intersubband mixing contribu- 
tion, then the above condition no longer holds. For this 
case we find (within a two-subband model) a critical value 
for the spin-orbit coupling (see Appendix A), which is 
given by fl so « (ira/W) 2 /[{it a/ W) + ak F ] — (3 c so , where 
Pso = t so /\t\, and k F is the Fermi wave number. The 
critical value [3 so will therefore define a weak spin-orbit 
coupling regime whenever /3 so < fi c so , and strong coupling 
regime if [3 so > fi c so . It should be noted that the critical 
value a c = ir 2 U 2 /[m*W(ir + Wk F )] of the Rashba spin 
orbit coupling parameter a that corresponds to f3^ a and 
therefore separates the weak and strong coupling regimes 
depends on the width W of the wire. 

In all our calculations, unless otherwise stated, the 
Rashba spin-orbit interaction is turned on and off adia- 
batically through a cosine-like function, with £ being the 
length of the adiabatic region in lattice sites. For the re- 
mainder of the paper we will work in units of \t\ for all the 
energies, hence we will refer to (3 so and t so interchange- 
ably. It is remarked that the calculations presented here 
are for narrow wires, in which the quantization in the 
transverse direction to the propagation is crucial. 

In Fig. 3 we show the ballistic conductance for four 
strength values of the spin-orbit hopping parameter t so 
for a narrow wire of width W = 60 nm and length 
L = 150 nm. Here as in the rest of the numerical re- 
sults and as discussed in Sec. II, purely spin-up polarized 
electrons are injected into the wire. Only the first mode 
is shown for clarity. In Fig. 3(a) the case t so = is plot- 
ted and since there is no electron spin-precession, all the 
transmitted electrons are spin-up polarized. When the 
spin-orbit parameter is increased to 0.03, Fig. 3(b), the 
precession effects become evident, about the 70% of the 
detected conductance is due to spin-up electrons, while 
about 30% is due to spin down electrons. For t so = 0.08, 
Fig. 3(c), the spin conductance is now reversed with re- 
spect to (a), that is, the net detected conductance is due 
only to spin-down electrons, = and G^ = 1.0 (in 
units of e 2 /h). Increasing the spin-orbit parameter fur- 
ther to 0.1, the spin- up conductance G^ ^ once again 
and the spin-down conductance G^ < 1.0, Fig. 3(d). 
Notice that in all cases, the spin conductance is almost 
independent of the incident Fermi energy; the situation 
is more complex when strong coupling is assumed, as 
we shall see later. The qualitative behavior of the spin 



G 



conductance described above is consistent (to a good ap- 
proximation) with the energy-independent electron spin 
modulation predicted by Datta and Das.El 

The spin conductance modulation is seen clearly when 
we plot GTJ. against the spin-orbit hopping parameter t so , 
see Fig. 4(a). Here the incident Fermi energy was fixed 
to 0.5 (kp w 0.7a -1 ) and W = 6a = 60 nm, which gives a 
critical value SO = 0.22. This value of j3 so separates the 
sinusoidal behavior of (predicted by Eqs. (13), (14) 
and (16)) for (3 so < 0.22 from its behavior for f3 so > 0.22 
(the strong spin subband mixing regime) where the con- 
finement energy is of the order of the intersubband mix- 
ing energy. The effect is clearer for a wire with W — 120 
nm, [see Fig. 4(c)] for which the critical value of (3 so is 
0.07. 

The importance of the intersubband mixing contribu- 
tion can be illustrated as follows. Let us write the Rashba 

(x) (y) 

Hamiltonian as the sum of two terms, H so = H so +Hso , 



with Hso^ = iu x <jyd / dx and Hso' = —ia y a x d/dy, where 
in general a x = a y = a for a 2DEG. Therefore, Eqs. 



Ay) 



(7) and (8) would be rewritten as {H so )}} m — 



a x (n\-^\m), and (H so )^ m = a y k y 5 nm - a x (n\-^\m), re- 
spectively. Now, since the only matrix elements that in- 
corporate the mixing of the different subbands are given 

by (n,a\Hso^\m,a'}, with n ^ m and a ^ a 1 , therefore 
by setting a x to zero (with a y ^ 0), the contribution 
of these matrix elements can be fully suppressed. This 
situation is shown in Fig. 4(b) and 4(d). Note that 
are very similar in Fig. 4(a) and 4(b) for t so < 0.22 
(t so < 0.07, case 4(c)) but not for t so > 0.22 (t so > 0.07, 
case 4(c)). Although setting a x — appears to be rather 
unphysical, it shows very explicitly, that the deviation 
from the sinusoidal modulation of G'^ for t so > 0.07 is 
owed essentially, to intersubband mixing induced by the 
Rashba spin-orbit coupling. 

The subband dispersion for t so — 0.08 is depicted in 
Fig. 5(a). The dispersion is calculated using the proce- 
dure described in the Appendix B, the rest of the sim- 
ulation parameters are the same as in Fig. 3. For this 
case, a linearly Rashba spin-split subband is obtained as 
expected , giving for this particular case Afc = 0.157 a -1 . 
The phase shift is then AO = AkL = it, where we have 
used the effective length of L — 20 a. Using now formula 
(16), we obtain for j3 so = 0.08 a phase shift AO = 3.2, 
just slightly above of what we obtain from the subband 
dispersion calculation£a The corresponding spin conduc- 
tance is dotted in Fig. 5(b). Note that a fully reversed 
spin-polarization occurs for all modes. Thus the spins 
precess in unison in all of the subbands and even a rela- 
tively wide multimode quantum wire should be expected 
to function as an efficient Datta-Das spin-transistor in 
this regime. For the same value of t so , Fig. 5(c) shows 
the conductance vs. the effective length L where the 
spin-orbit coupling is present. The oscillations of G^ and 
have a differential phase which corresponds exactly to 
the expected value for length L = 20 a in Eq. (16). It is 
important to emphasize that the qualitative features seen 



here in the conductance for weak t so remain the same for 
wider wires and with a lateral constriction in the region 
with spin-orbit coupling. For example in Fig. 6(a) we 
have plotted the conductance for a wire having twice the 
width discussed above, i.e. here W — 120 nm, where 
t so = 0.05 and with a strength for a parabolic confining 
potential of w x — 0.2 (here the confinement potential is 
given by w x (x/a) 2 , with w x in units of |t|) . We note 
the same step-like characteristic of the ballistic conduc- 
tance. Note that here as well, (Fig 6(a))the modulation 
achieved between G^ and G^ remains constant regardless 
of the incident Fermi energy and number of populated 
subbands. 

Now we turn our attention to the case [3 so > 
(ira/W) 2 /[(na/W) + akp], i.e. strong spin-orbit inter- 
action, for the case with W = 60 nm, (3 so > 0.22. Fig. 
7 displays the spin conductance for four different val- 
ues of the strength for t so - Here the adiabatic region 
is set to £ = 20a, and the rest of the parameters are 
as in Fig. 3. For sake of simplicity we focus only on the 
first propagating mode. It is evident that for this regime, 
the spin-conductance behavior is markedly different from 
that with weak spin-orbit coupling discussed previously. 
Here G^i are both strongly energy dependent. We note 
that the probability of an injected spin up electron to 
pass unchanged through the spin-orbit region decreases 
with energy [Fig. 7(a)], while the probability of detecting 
a transmitted spin down electron increases accordingly. 
For t so = 0.3, Fig. 7(b), a full polarization develops 
when reaching an incident energy of Ep = 0.8, although 
without the spin flipping. Cases (c) and (d) are quite 
interesting. A clear precession of G'* with energy is ob- 
served, having a sinusoidal like character. For instance, 
in (c), the polarization of the transmitted current has 
opposite orientation than the polarization of the injected 
current for Ep = 0.65, and is reversed again at Ep = 1.0. 
This surprising result (which is due to intersubband cou- 
pling and is therefore beyond the scope of the original 
theory of Datta and Das) would suggest that for certain 
values of the strength of the spin-orbit interaction, just 
by varying the the Fermi energy, the device can work as 
a spin "switch" device, while the spin-orbit interaction is 
kept constant. In other words, it can function as a spin- 
transistor such that the switching can be tuned just by 
varying the Fermi energy rather than varying the Rashba 
coupling, as in Datta and Das spin-transistor. 

The subband dispersion for the case t so = 0.4 of Fig. 
6(b) is shown in Fig. 8(a). For this large value of t so , the 
spin-subband dispersion deviates considerably from the 
typical linear Rashba splitting for the weak spin-orbit 
interaction (A similar-dispersion was obtained recently 
by Barnes and Morozo). The subbands are no longer 
parabolic which gives rise to a dependence of Afc on the 
energy. Here Afc is the difference between the wave vec- 
tors associated the two spin subbands at a given energy, 
Afc -> Ak{E) = k 2 (E) - h(E). We find for example 
that Afc(£ F = 0.5) = 0.6489 a" 1 , whereas Afc(£ F 
I .() = 0.4533 a" 1 , which yields A6(E F = 0.5) = 6.2?r 
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and A9(Ep — 1.0) = 4.37T, where we have used an effec- 
tive length L — 30 a. The net change in the differential 
phase between these two cases is thus 5(A9) « 2ir. This 
last result can be independently checked by looking at 
the spin-resolved conductances G^ and in Fig. 8(b) 
where the oscillations of both G^ and G^- exhibit phase 
shifts of ~ 2ir between Ep — 0.5 and Ep = 1.0. We have 
also plotted the spin conductance for Ep = 1.0 versus 
the length of the spin-orbit interaction region in the wire, 
[Fig. 8(c)]. A wave length A of 13.8 a, is extracted from 
the data. The associated wave number is Afc = 0.45, 
value that matches that obtained independently from the 
band dispersion calculation. 

Now let us return to Fig. 6 in order to analyze the 
strong t so regime for a wider wire (L = 12a). In partic- 
ular Fig. 6(b) shows an analogous case to that studied 
in Fig. 6(a) but with t so = 0.4. It is evident that the 
strong coupling here suppresses the plateaus for Gfi+ to 
a significant degree. A rather complicated structure is 
obtained revealing the non-trivial nature of the subband 
intermixing. We have also calculated the relative conduc- 
tance change [defined as AG/G Q = (G T -G l )/(G^ +G 1 )} 
against the spin-orbit parameter t so for three different 
effective lengths of the wire with the Rashba interaction 
(L = 4W, 8W, and 15W). This is shown in Fig. 9 for 
a wire with width W = 6a and uniform parabolic po- 
tential strength of w x — 1.8 over an effective length of 
10a and at the Fermi energy Ep — 2.5 for the three 
cases. A beat-like pattern is found due to the spin pre- 
cession, developing nodes as the length of the channel is 
increased.LJ In the inset the channel length dependence 
of AG/G is plotted for a typical experimental t so , it 
is clear that for such relatively weak spin-orbit interac- 
tion the relative conductance change has a negative slope 
with a change in sign (indicative of a spin-precession) in 
the length range shown. This behavior resembles that 



observed recently by C. M. Hu et alE3 in their length 
dependence measurements for the spin precession. 

A final comment on the adiabaticity: the adiabatic- 
ity in the spin-orbit interaction was introduced in our 
calculations in order to model a smooth transition be- 
tween the regions with no spin-orbit coupling (i.e near 
the ferromagnetic contacts) and the region with the fi- 
nite Rashba spin-orbit coupling along the quantum wire. 
For the weak coupling regime the calculations of the 
spin-conductance vs. Fermi energy showed plateaus with 
small oscillations in the non-adiabatic case, nevertheless 
the spin-precession behavior was found to be qualita- 
tively very similar to that observed in the adiabatic case. 
However, the strong coupling regime showed drastic dif- 
ferences. Calculations without the adiabaticity showed a 
rather complicated behavior with a significant suppres- 
sion of the conductance plateaus. 

In summary, we have shown that a strong Rashba spin- 
orbit interaction can produce dramatic changes in the 
spin-resolved transmission of spin-polarized electrons in- 
jected into ballistic narrow wires. The effects can be 
very significant, and can even suppress the expected spin- 
modulation, as the strong (Rashba) coupling was found 
to induce a complex dependence of the spin-precession 
on the incident Fermi energy of the injected electrons. 
These results should be of importance for the spin- 
injection (modulation) in quasi-one-dimensional devices 
under large bias used to tune the Rashba effect. 
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APPENDIX A: CRITERIA FOR WEAK AND STRONG RASHBA COUPLING 



In this appendix the criteria we have used to distinguish the regimes for weak and strong Rashba spin-orbit coupling 
in a Q1DEG system are deduced analytically. Consider |<j>) the eigenstate solution of 



(Al) 



where H = H a + H so is the effective mass Hamiltonian for a Q1DEG in the presence of the Rashba coupling, with 



H \n,a) = E%\n,a), and E° = E n + h z fc2/2m*, ( n is the subband index and a =|, j), see Sec. II. B. We expand |$) 
in terms of the eigenstates for the system in the absence of the Rashba effect. That is |$) = J2 n a a n<j\ n , o - }- Inserting 
this expansion into (Al) and multiplying from the left by {m,a'\ yields, 



(A2) 



Now since (H so ) 



a. a 

nm 



for all a = a' (sec cq. (9)), therefore eq. (A2) will read 

( E n - E ) a n<7 + ( H so)nn a na> + (4o)nm W = ^, 



(A3) 



the second term in (A3) which is diagonal in the subband index will corresponds to the linear Rashba splitting, while 
the third term gives rise to the spin subband intermixing. In a two-subband model (denoted below by n and n+1), 
the system of eqs. (A3) can be rewritten as 



(E° n -E ak y 



V 



ak y E° - E 

aA so 
-aA so 




(A4) 



with the definition of the intersubband mixing term, A so = (n + l\4-\n). The eigenvalues of (A4) are explicitly. 



E 1 ,2 = -(E° n+1 +E° n )±-^/(AE w 



2ak y ) 2 + (2aA so ) 2 , 



and 



E 3A - l -{E° n+1 + E° n ) ± ^(AE w + 2ak y r + (2aA s 



(A5) 



(A6) 



with AEw = E° +1 — E° = E n+ \ — E n . Notice that for small aA so , after a zeroth order Taylor-series expansion of 
the radicals in (A5) and (A6), the eqs. will reduce exactly to that given by eq. (12) for the energy subband linear 
splitting. Such expansion will hold as long the condition 



(2aA s 



(AE W ± 2ak v ) 2 



< 1. 



(A7) 



is satisfied. Thus, Eq. (A7) will lead us to a criteria for which the subband intermixing is, on the other hand, not 
neglectable. For such case, explicitly 



2aA s 



AEw ± 2ak y 



1. 



(A8) 



Modeling the potential that defines the Q1DEG as a hard well confining potential, we have AEw = E° +1 — E^ = 

(2n + l)(fpF)|t|, and A so = "^"ff (w 7 )- Inserting this expressions in (A8) and replacing ky — ► hp, where kp is the 
Fermi wave number, we get 



a 



(2n + l) 2 (na) 2 \t\ 



8n(n + l)W ± 2(2n + l)k F W 2 



(A9) 



therefore, for the first subband (n = 1), and using a — 2af3 so \t\, we get (choosing the + solution in order to get the 
lowest value for the ratio in (A9)) 
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{ira/W) 2 
(ira/W) + aki 



(AlO) 



Eq. (AlO) defines the critical value of the spin-orbit coupling at which the intersubband mixing becomes relevant. 
Therefore we can define a regime where the spin-orbit coupling is weak, whenever f3 so < (ira/W) 2 /[(wa/W) + akp], 
and a regime where the spin-orbit coupling is strong, such that spin subband intermixing becomes important, that is 
whenever (3 so > (na/W) 2 /[(na/W) + ak F ]. 



APPENDIX B: SUBBAND SPECTRUM CALCULATION 



In this appendix we describe the transfer matrix method for the calculation of the subband structure of an uniform 
quantum wire with a Rashba spin-orbit interaction. The method is a generalization of the spin-independent Ando 
calculation.^ The dispersion is calculated within the lattice model for a wire with a finite and uniform spin-orbit 
interaction and a lateral confinement potential. 

By applying the states (lma\ and Y$>) to the left and right of Eq. (19), respectively, we arrive at the following 
discrete coupled equations [for clarity, the notation a — ±, is used instead of the arrow symbols to denote spin-up 
and spin-down states] 



( W tm ~ E )xtrn + 1 (xf +1 , m + xt-l,m + X^m+1 + xtn-l) T *~ (xf+1,- 



,m 



it 



Xl m — 1 



= 
(Bl) 



Here we defined Wf m = e + Vi. m . Now, if we set the on-site energy e = — it, and since a uniform (along the y-axis) 
confinement potential is assumed, Vi, m — > Vi, and thus, W° m — > Wi. Hence the equation above can be written in a 
matrix form as follows 



(H G E)A^ + T^ A^ + t(A, m+1 + A m _ 1 ) + it so (A^ n _ 1 A^ +1 ) — 



Here we have defined the vector A!L as the rath column 



(B2) 



/ Xl,m \ 
X2,m 

Xl t m 

\ XN^,m J 



(B3) 



and the H G and T^ are N x x N x matrices given by 
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(B4) 



Now, since both the confinement potential and the spin-orbit interaction are assumed to be uniform along the 
longitudinal direction (propagation direction), the system has translational symmetry along the y-axis. Therefore we 
can make use of the Bloch theorem. We can define then A^ +1 = e lkya A^ n and A a m — e lkya A^ n _ 1 , hence A^ +1 = AA^, 
with A = e lky °". Inserting this last expression for Af n+1 in Eq. (A2), we get 



(H - E)At + T±A* + AiA± 
which can be rewritten in a matrix form as follows, 



tA 



m — 1 



(B5) 
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/ H 



V 



o-E T+ 


t 


it so 


T so H o - E 


it 
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1 








1 









/ A 4 
A+ 



= A 



-t 












-t 














1 














1 



/ A 4 

K 

A+ 



(B6) 



■1/ 



By solving the (4A^) x (4N X ) generalized eigenvalue problem, for a given incident Fermi energy E, all the eigenvalues 
A are determined, and hence all k y (E), yielding the desired subband spectrum. 
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FIG. 1. (a) Cross sectional schematic of Datta-Das spin modulator device, (b) cross sectional schematic of a split gate device 
used to create a Q1DEG, (c) diagram of the quasi-one dimensional quantum channel. 
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FIG. 2. (a) Schematic of the tight-binding model for the system. In the shaded areas the spin-orbit interaction is finite, 
tso 7^ 0. (b) Schematic of a quasi-one dimensional spin modulator device. 
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FIG. 3. Ballistic spin-conductance in units of e 2 /h against the Fermi energy for a narrow quantum wire of width W = 6a, 
showing the spin-modulation for different values of the strength for the spin-orbit hopping parameter. Filled circles G T , empty 
circles G^, and solid line G T + G l . In all cases (as in the rest of the calculations), exclusively spin-up polarized electrons are 
injected to the wire. Note that the spin-conductance modulation is almost independent of the incident Fermi energy as in the 
model of Datta and Das (Ref. [1]). 
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FIG. 4. Spin-orbit coupling strength dependence of the ballistic spin-conductance, solid line G T , dashed line G^: (a) Narrow 
wire of W = 6a, and uniform spin-orbit coupling (a x = a y — 2at so ). (b) Same as in (a) but with a x = and a y — 2at so , 
perfect oscillations are seen for all t so . (c) Same as in (a) with W = 12a. (d) Modulation for W = 12a, with a x = and 
a y = 2at so . The intersubband mixing clearly changes the otherwise perfectly sinusoidal spin-conductance modulation. 
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FIG. 5. (a) Energy subband dispersion for a narrow wire of W = 6a and t so = 0.08 (weak spin-orbit coupling regime), 
(b) Ballistic spin-conductance for the wire in (a), note that for this value of t so , practically all the conductance is due to the 
precessed spin-down electrons (dashed line), while that for the spin-up is almost zero for all energies (solid line), (c) Oscillating 
behavior of the spin-conductance versus the effective length of the spin-orbit interaction region, here t so = 0.08 as well. A wave 
length of A = 20a is extracted from the data. 
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FIG. 6. Fermi energy dependence of the ballistic spin-conductance (in units of e 2 /h) for a wire of width W 
confinement strength (curvature) of the parabolic potential w x = 0.2. (a) Weak spin-orbit coupling case, t so — 0.05. 
spin-orbit coupling case, t so = 0.4. Filled circles G T , empty circles , and solid line G T + G^. 
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FIG. 7. Ballistic spin-conductance for the strong spin-orbit coupling regime, (t ao > 0.22) for the wire of Fig. 3. Only the 
first propagating mode for each spin is shown. In this regime, both G T and G l are strongly dependent on the incident Fermi 
energy, with an oscillating behavior as t so is enhanced, which exhibits the energy dependence of the spin-precession. Cases (c) 
and (d) are particularly interesting. For Ef ~ 1.0 in (d) the polarity of the transmitted electrons is reversed (with respect 
to the injected spin-up polarized electrons) Perfectly polarized spin up or spin down electrons can be emitted from the wire 
depending on the Fermi energy. 




FIG. 8. (a) Energy subband dispersion for a narrow wire of W — 6a and t so = 0.4 (strong spin-orbit coupling regime). Note 
that the subbands are no longer parabolic as in weak t so regime, resulting in a Ak — » Ak(E). (b) Ballistic spin-conductance for 
the wire in (a); filled circles G\ empty circles G , and solid line G T + G . The energy dependence is clearly evident here, (c) 
Oscillating behavior of the spin-conductance versus the effective length of the spin-orbit interaction region (t so = 0.4) . The 
wave length of the oscillations is A ?a 13. 8a. 
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FIG. 9. Plots of the relative conductance change [AG/G = (G T — G- L )/(G T + G 1 )] against the strength of spin-orbit parameter 
for three different effective lengths of the spin-orbit (Rashba) interaction, L — 4W, L = SW, and L = 15W, with W = 6a. 
Clearly the beat-like behavior is enhanced with L, showing the nature of the electron spin-precession induced by the Rashba 
effect. In the inset of (a) the relative conductance change is plotted as a function of the effective length for a weak spin-orbit 
hopping parameter, t so — 0.015. A change in sign for AG/G as the length is increased is evident showing that the effect of 
the spin-precession still important for such regime. 
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